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ABSTRACT 

In this paper, we study the downlink multiuser scheduling 
problem for systems with simultaneous wireless information 
and power transfer (SWIFT). We design optimal scheduling 
algorithms that maximize the long-term average system 
throughput under different fairness requirements, such as 
proportional fairness and equal throughput fairness. In par¬ 
ticular, the algorithm designs are formulated as non-convex 
optimization problems which take into account the minimum 
required average sum harvested energy in the system. The 
problems are solved by using convex optimization techniques 
and the proposed optimization framework reveals the tradeoff 
between the long-term average system throughput and the 
sum harvested energy in multiuser systems with fairness 
constraints. Simulation results demonstrate that substantial 
performance gains can be achieved by the proposed optimiza¬ 
tion framework compared to existing suboptimal scheduling 
algorithms from the literature. 

Index Terms — RF energy harvesting, wireless informa¬ 
tion and power transfer, optimal multiuser scheduling. 

1. INTRODUCTION 

Over the past decades, battery-powered devices have been 
deployed in many wireless communication networks. How¬ 
ever, since batteries have limited energy storage capacity and 
their replacement can be costly or even infeasible, harvesting 
energy from the environment provides a viable solution for 
prolonging the network lifetime. Although conventional 
natural energy resources, such as solar and wind energy, 
are perpetual, they are weather-dependent and location- 
dependent, which may not suitable for mobile communication 
devices. Alternatively, background radio frequency (RF) 
signals from ambient transmitters are also an abundant source 
of energy for energy harvesting (EH). Unlike the natural 
energy sources, RF energy is weather-independent and can 
be available on demand. Nowadays, EH circuits are able 
to harvest microwatt to milliwatt of power over the range 
of several meters for a transmit power of 1 Watt and a 
carrier frequency less than 1 GHz Q. Thus, RE energy 
can be a viable energy source for devices with low-power 
consumption, e.g. wireless sensors @ 0 - Moreover, RE EH 
provides the possibility for simultaneous wireless information 
and power transfer (SWIFT) since RF signals carry both 
information and energy 0§. 

The integration of RF EH into communication systems 
introduces a paradigm shift in system and resource alio- 
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cation algorithm design. A fundamental tradeoff between 
information and energy transfer rates was studied in @16). 
However, current practical RF EH circuits are not yet able 
to harvest energy from an RF signal which was already used 
for information decoding (ID) @. To facilitate simultaneous 
ID and EH, a power splitting receiver was proposed in 0 
and @. The energy efficiency of a communication system 
with power splitting receivers was investigated in @|. In 
addition, a simple time-switching receiver has been proposed 
which switches between ID and EH in time. Furthermore, 
multiuser multiple input single output SWIFT systems were 
studied in 1T0|, where beamformers were optimized for 
maximization of the sum harvested energy under minimum 
required signal-to-interference-plus-noise ratio constraints 
for multiple ID receivers. In 1ID> the optimal energy 
transfer downlink duration was optimized to maximize the 
uplink average information transmission rate. In GD- 
GD, beamforming design was studied for secure SWIFT 
networks with different system configurations. In p6) , a 
multiuser time-division-multiple-access system with energy 
transfer in the downlink (DL) and information transfer in 
the uplink was studied. The authors proposed a protocol 
for sum-throughput maximization and enhanced it by fair 
rate allocation among users with different channel conditions. 
Nevertheless, multiuser scheduling, which exploits multiuser 
diversity for improving the system performance of multiuser 
systems, has not been considered in 0-G3- Recently, 
simple suboptimal order-based schemes were proposed to 
balance the tradeoff between the users’ ergodic achievable 
rates and their average amounts of harvested energy in fTT) . 
However, the scheduling schemes proposed in | fT7| are unable 
to guarantee quality of service with respect to the minimum 
energy transfer. In fact, optimal multiuser scheduling 
schemes that guarantee a long-term minimum harvested 
energy for SWIFT systems have not been considered in the 
literature so far. 

Motivated by the above observations, we study optimal 
scheduling schemes for long-term optimization which control 
the rate-energy (R-E) tradeoff under the consideration of 
proportional fairness and equal throughput fairness. 

2. SYSTEM MODEL 

We consider a SWIFT system that consists of one access point 
(AF) with a fixed power supply and N battery-powered user 
terminals (UTs), see Fig.@ The AF and the UTs are equipped 
with single antennas. Besides, we adopt time-switching 
receivers at the UTs @ to ensure low hardware complexity. 
We study the user scheduling for DL transmission. We 
assume that the transmission is divided into T time slots and 
in each time slot perfect CSI is available at the AF. Also, the 
data buffer for the users at the AF is always full such that 
enough data packets are available for transmission for every 



where 



Fig. 1. A multiuser system with SWIFT for TV = 8 time¬ 
switching user terminals (UTs). 

scheduled UT. In each time slot, the AP schedules one user 
for ID, while the remaining users opportunistically harvest 
energjQ from the received signal. We assume block fading 
channels. In particular, the channels remain constant during a 
time slot and change independently over different time slots. 
Besides, the users are physically separated from one another 
such that they experience independent fading. 

Furthermore, we adopt the EH receiver model from |j7). 
The RF energy harvested by user n S {1,..., N} in time slot 
iG{l,...,r}is given by 

Qn{i) = ^nPhn{i), (1) 

where P is the constant AP transmit power, 0 < Cn < 1 
is the RF-to-direct-current (DC) conversion efficiency of the 
EH receiver of user n, and Ti„(z) is the channel power gain 
between the AP and user n in time slot i. 

3. OPTIMAL MULTIUSER SCHEDULING 

In the following, we propose three optimal multiuser schedul¬ 
ing schemes that control the R-E tradeoff under different 
fairness considerations. 

3.1. Optimal Maximum Throughput (MT) Scheme 

First, we consider a scheduling scheme which maximizes the 
average sum rate subject to a constraint on the minimum 
required average aggregate harvested energy. We note that 
this scheme aims to reveal the best system performance, and 
fairness in resource allocation for UTs is not considered. 
To facilitate the following presentation, we introduce the 
user selection variables q„(i), where i G {1,2, ...T} and 
n G {!,..., N}. In time slot i, if user n is scheduled to 
perform ID, q„(i) = 1, whereas = 0, Vn ^ n, i.e., all 
the remaining idle users harvest energy from the transmitted 
signal. Now, we formulate the MT optimization problem as 
follows. 

Problem 1. Maximum Throughput Optimization: 
maximize Rsum 

qn{i) 

N 

subject to Cl; ^g„(i) = l,Vt, ^2) 

n—1 

C2: € {0, l},Vn,z, 

C3. (^sum ^ ^req: 

*We consider a unit-length time slot, hence the terms “power” and 
“energy” can be used interchangeably. 


^ T N 

Rsum — ^ ^ ^ ^ (3) 

T^oo 1 

n—1 

^ T N 

Osum = lim y^(l - qn{i))Qn{i), and (4) 

T^oo 1 

i—1 n—1 

Cnii) = l0g2 ^ . (5) 

Here, cr^ is the additive white Gaussian noise power at UT 
n. In the considered problem, we focus on the long-term 
system performance for T ^ oo. Constraints Cl and C2 
ensure that in each time slot only one user is selected to 
receive information. C3 ensures that the average amount of 
harvested energy Qsum is no less than the minimum required 
amount Qreq- Since the user selection variables Vn, *, 
are binary, problem (|^ is non-convex. In order to handle 
the non-convexity, we adopt the time-sharing relaxation. In 
particular, we relax the binary constraint C2 such that qn{i) is 
a continuous value between zero and one. Then, the relaxed 
version of problem Q can be written in minimization form 
as: 

minimize — i?sum 

gri(i),Vn,2 

subject to Cl, C3, (6) 

C2 : 0 < qn{i) < l,Vn,T 

Now, we introduce the following theorem that reveals the 
tightness of the binary constraint relaxation. 

Theorem 1. Problems ® and (|^ are equivalent with 
probability one, when /i„(I),Vn,i are independent and con¬ 
tinuously distributed. In particular, the constraint relaxation 
of C2 is tight, i.e., 

C2 C2 : 0 < < 1, Vn, z. (7) 

Proof. Theorem 1 will be proved in the following based on 
the optimal solution of •S- □ 

In other words, we can solve 0 via solving It can 
be verified that the relaxed problem is convex with respect 
to the relaxed optimization variables and satisfies the Slater’s 
constraint qualification. Therefore, strong duality holds and 
the optimal solution of (|^ is equal to the optimal solution 
of its dual problem ID- Thus, we solve Hi via the dual 
problem. To this end, we first define the Lagrangian function 
for the above optimization problem as 

P{qn (T), (0 ; 7n (0 ; -^sum 

T / N \ T N 

+1 ) + XIXI 

2—1 \n —1 / 2—1 n —1 

T N 

-XX/^ n(^)Q'n(^) ^ (Qreq Q sum) ; (8) 

2=1 n=l 

where A(i),{a„(z) > 0,(3n{i) > 0}, and z/ > 0 are the 
Lagrange multipliers corresponding to constraints Cl, C2, 
and C3, respectively. Thus, the dual problem of (j^ is given 
by 

maximize minimize A(z), a„(z),/3„(z), z/). (9) 

«r.(j).4n(*)>0,A(i) q„{i) 

^Here, “equivalent” means that both problems share the same optimal 







In order to determine the optimal user selection policy, we ap¬ 
ply standard convex optimization techniques and the Karush- 
Kuhn-Tucker (KKT) conditions. Thereby, we differentiate 
the Lagrangian in ([^ with respect to g„(i) and set it to zero 
which yields: 

dL{...) Qnii) Cn{i) . 

a = H^)+au{^)-|3n{^)+v^f; -^ = 0, Vn,z. 

OQnW T T 

( 10 ) 

We define n* as the optimal user selection index for ID at 
time slot i, i.e., = 1 and qn{i) = 0, Vn ^ n*. From 

the complementary slackness condition, we obtain Q!„(z) = 
0, Vn ^ n* and /3„. (z) = 0. Now, we denote the optimal 
dual variable for constraint C3 as v* and substitute it into ( [T0| i. 
Then, the selection metric for UT n is given as 

= T(A(z)+a„(i) - /3„(i)) = C„(f)-zz*Q„(i). 

( ll a) 

Hence, the selection metric for the scheduled UT is 
A“T(z) = T (A(z) + (z)) = (z) - (i), 

( ll b) 

and the selection metric for the remaining UTs is 

A“T(*) = T (A(z) - /3„(z)) = C„(z) - u*Qn{i), Vn^n*. 

( ll c) 

Subtracting (|1 lc|i from (|1 lb|i yields 


Problem 2. Optimal Proportional Fair Optimization: 
minimize — i?sum 

qn{i)yn^i 

subject to Cl, C2, C3, 

T 

C4: ^ = 0,yn, 

i=l 

where C4 specihes that each UT has to access the channel 
for ^ number of time slots. For the tightness of the binary 
relaxation, please refer to Theorem 1. 

Now, we solve ( |T4| ) via convex optimization techniques 
by following a similar approach as in the previous section. 
The Lagrangian function for problem ( [l4| l is given by 

(^) 5 ^(^); (^) 7 /^n (f) ; — .(^suni 

T / N \ T N 

+ 1 ) + XIXI 

2—1 \n —1 / 2—1 n —1 

T N 

-XX/^ n(^)9n(^) “I” ^ (Qreq Q sum) 

2^1 n —1 

^ ^ 1 \ 

+ X^" 

n=l \ 2=1 / 


A“T(z) - Ar (*) = + /3n{^)) ■ (12) 

Since a„-(z) + /?«(*) > 0 from the dual feasibility 

conditions, we obtain Aj^J'(z) > Ajf"'"(z), Vn ^ n*. 
Furthermore, Ajf^(z) Vn, are continuous random variables, 
therefore Pr{Aj^.'^(z) = Ajf"'"(z)} = 0, Vz, where Pr{-} 
denotes the probability of an event. Thus, A^.^(z) > 

A“T(i), Vn ^ n* and the optimal selection criterion for the 
MT scheme in time slot i reduces to 

f 1 ifA!f"'’(z)= max Aj^'^(z) 

qn[i) = { ‘ 'V (13) 

0 otherwise 


In other words, the solution of the relaxed problem is itself of 
the Boolean type. Therefore, the adopted binary relaxation 
is tight. Besides, v* depends only on the statistics of the 
channels. Hence, it can be calculated offline, e.g. using the 
gradient method, and then used for online scheduling as long 
as the channel statistics remain unchanged. We emphasize 
that although the original problem in (|^ considers inhnite 
number of time slots and long-term averages for the sum rate 
and the total harvested energy, the optimal scheduling rule 
in ( [TJI depends only on the current time slot, i.e., online 
scheduling is optimal. 


3.2. Optimal Proportional Fair (PF) Scheme 

In the MT scheme, UTs with weak channel conditions may 
be deprived from gaining access to the channel which leads 
to user starvation. In order to strike a balance between system 
throughput and fairness, we introduce proportional fairness 
into our scheduler, which aims to provide each UT with a 
performance proportional to its channel conditions. This is 
achieved by allowing all UTs to access the channel with equal 
chances. In this case, the optimization problem with the 
relaxed binary constraint on the user selection variables is 
formulated as: 


where A(z),{a„(z) > 0,/3„(z) > 0},z^ > 0, and 7 „ are 
the Lagrange multipliers corresponding to constraints Cl, C2, 
C3, and C4, respectively. By using the KKT conditions, we 
obtain the following UT selection metric: 

APP(z) = C„(z) - iz*Q„(z) - 7 ^, (16) 


where the optimal Lagrange multipliers 7 * ensure that each 
user accesses the channel on average an equal number of 
times. Thus, the optimal selection criterion for the PF scheme 
is 


qn{i) 


1 if API" (i) =^max^^ 
0 otherwise 


(17) 


We note that the optimal PF scheduling rule is similar to the 
MT scheduling rule in ( [T3] l, but the PF seelction metric in 
(fThll contains an additional term 7 * that provides proportional 
fairness. Also, v* and 7 * can be calculated offline using the 
gradient method. 


3.3. Optimal Equal Throughput (ET) Scheme 

Although the PF scheduler enables equal channel access 
probability for all UTs, it does not provide any guaranteed 
minimum data rate to them. On the contrary, the ET criterion 
is more fair from the users’ prospective compared to the PF 
criterion, as all the UTs achieve the same average throughput 
asymptotically for T —>■ 00 . Therefore, in this section, we 
design a scheduler which achieves ET fairness. Thus, the 
objective is to maximize the minimum average achievable 
rates among all the UTs, i.e., maximize minC^ where 

n 

Cn = limT->oo ^ qn{i)Cn{i)- Using Theorem 1 , we 
formulate our equivalent convex optimization problem in its 
hypograph form. 








Problem 3. Optimal Equal Throughput Optimization: 
minimize — r 

r,qji{i)yn,i 

subject to C1,C2, C3, 

C5; r - < 0, Vn, 

where r is an auxiliary variable. The Lagrangian function for 
problem in ( [T8] l is given by 

L{qnii),X{i),an{i),/3nii),i',0n) = -r 

T / N \ T N 

i—l \n—l / i—1 n—1 

T N N 

-EE/3 n{i)qn{i)+l' (Qreq- Qsum) E 9n (r - C„) , 

i—l n—1 n—1 

(19) 

where A(i), {a„(j) > 0, /3„(i) > 0}, > 0, and > 0 are 

the Lagrange multipliers corresponding to constraints Cl, C2, 
C3, and C5, respectively. By using the KKT conditions, we 
obtain the UT selection metric for ET scheduling: 

A^^(i)=CC„(z)-iy*Q„(z), (20) 

where the optimal Lagrange multipliers 0* ensure that all 
users have ET. Thus, the optimal selection criterion for the 
ET scheme is given by 

1 ifA®'^(i)= max 

‘ . ( 21 ) 

0 otherwise 

Again, the gradient method can be used to obtain the optimal 
values for ly* and 0* offline by utilizing the channel statistics. 

Remark 1. We note that the above considered problems can 
be formulated as Markov Decision Process (MPD) or solved 
via Lyapunov optimization approach, please refer to l{19^ for 
details. 

4. SIMULATION RESULTS 

In this section, we evaluate the performance of the proposed 
scheduling schemes using simulations. The important sim¬ 
ulation parameters are summarized in Table 1. We adopt 
the path loss model from | |20| and the UTs are randomly 
and uniformly distributed between the reference distance and 
maximum service distance. Eor comparison, we also show the 
performance of the following suboptimal scheduling schemes 
from i fTTl : 

1. Order-based MT scheduler: The scheduling rule is 
n*(i) = argorder hn{i), where argorder is defined as 

nG{l....,W} 

the argument of a certain selection order j S {1,..., N}. 
In other words, the user whose channel power gain hn{i) 
has order j is scheduled for ID. 

2. Order-based PE scheduler: The scheduling rule is n*{i) = 

argorder where fin denotes the mean channel 

power gain of UT n. 

3. Order-based ET scheduler: The scheduling rule is 
n*(i) = arg min rn(i — 1), where Ou„ G {1) • ■ • 

Ou„ GiSa 


Table 1. Simulation parameters. 


Parameter 

Value 

AP transmit power P 

40 dBm 

Noise power cr^ 

—62 dBm 

RP-to-DC conversion efficiency 

0.5 

Path loss exponent 

3.6 

Maximum service distance 

100 m 

Reference distance 

2 m 

Antenna gain of AP and UTs 

10 dBi & 2 dBi 

Carrier center frequency 

915 MHz 

Bandwidth 

200 kHz 

Eading channel 

Rayleigh 



Fig. 2. Average sum rate versus average sum harvested energy 
of the MT schemes for different numbers of UTs. 

is the order of the instantaneous normalized signal-to- 
noise-ratio of user n, 5a is a predefined set of orders, 
where only the users set Ou„ G { 1 ,..., A^} fall into 
5a are eligible for being scheduled, and r„(* — 1 ) is the 
throughput of user n averaged over all previous time slots 
up to time slot i — l. 

Fig. |2] shows the average sum rate (bits/(channel use)) 
versus the average sum harvested energy (Watts) of the MT 
schemes for different numbers of users. We note that the 
suboptimal order-based scheme can only achieve discrete 
points on the R-E curves, corresponding to the selection 
orders j G A^}. On the contrary, the proposed optimal 

MT scheduling scheme can achieve any feasible point on 
the R-E curve, which provides a higher flexibility for the 
system designer to strike a balance between average sum 
rate and average harvested energy. Besides, as expected, the 
average system sum rate increases with the number of UTs 
as the proposed scheme is able to exploit multiuser diversity. 
Eurthermore, the average sum harvested energy also increases 
with the number of UTs since more idle users participate in 
energy harvesting in any given time slot. 

Eig. 1^ and Eig. depict the average sum rate 
(bits/(channel use)) versus the average sum harvested energy 
(Watts) for the PE and ET schemes, respectively. It can be 
seen that the feasible R-E region of all schemes decreases 
compared to the MT scheduler in Eig. This is because 

















Fig. 3. Average sum rate versus average sum harvested energy Fig. 4. Average sum rate versus average sum harvested energy 

of the PF schemes for different numbers of UTs. of the ET schemes for different numbers of UTs. 


both the PF and the ET schedulers take fairness into account 
in the resource allocation and, as a result, cannot fully exploit 
the multiuser diversity for improving the average system sum 
rate. On the other hand, it can be seen that our proposed 
optimal schemes provide a substantial average sum rate 
gain compared to the corresponding suboptimal order-based 
schemes, especially for a high amount of average harvested 
energy in the system. In fact, the proposed optimization 
framework provides more degrees of freedom across different 
time slots in resource allocation compared to the suboptimal 
scheduling schemes. This allows the system to exploit the 
multiuser diversity to some extent for resource allocation even 
if fairness is taken into consideration. 

5. CONCLUSION 

In this paper, we have proposed optimal multiuser scheduling 
schemes for SWIPT systems considering different notions 
of fairness in resource allocation. The designed schemes 
enable the control of the tradeoff between the average sum 
rate and the average amount of sum harvested energy. Our 
results reveal that for the maximization of the system sum rate 
with or without fairness constraints, the optimal scheduling 
algorithm requires only causal instantaneous and statistical 
channel knowledge. Simulation results revealed that sub¬ 
stantial performance gains can be achieved by the proposed 
optimization framework compared to existing suboptimal 
scheduling schemes. 
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